Melting transition of an Ising glass driven by magnetic field. 
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The quantum critical behavior of the Ising glass in a magnetic field is investigated. We focus on 
the spin glass to paramagnet transition of the transverse degrees of freedom in the presence of finite 
longitudinal field. We use two complementary techniques, the Landau theory close to the T = Q 
transition and the exact diagonalization method for finite systems. This allows us to estimate the 
size of the critical region and characterize various crossover regimes. An unexpectedly small energy 
scale on the disordered side of the critical line is found, and its possible relevance to experiments on 
metallic glasses is briefiy discussed. 
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Understanding disordered systems is one of the main 
challenges of condensed matter physics, since the pres- 
ence of disorder is always unavoidable in experiments. 
When disorder is strong it can dominate the physics 
and lead to exotic states of matter such as the glassy 
phases 0- The most salient properties observed in glassy 
systems are the slow dynamical relaxation and history 
dependence of thermodynamics. Research on quantum 
spin systems is of primary importance because of poten- 
tial technological applications. Current work in quantum 
computing and spintronics, where the understanding of 
relaxation processes is crucial 0,0, is boosting a renewed 
interest in basic models of disordered quantum magnets. 

The goal of the present work is to consider the ran- 
dom Ising model that displays a quantum paramagnet 
to spin glass transition driven by fluctuations introduced 
by an external magnetic field. We tackle the problem by 
utilizing two different theoretical approaches. We solve 
the model using the recently introduced technique of ex- 
act diagonalization that includes the averaging over an 
ensemble of disorder realizations in a flnite system. The 
relevant results arc then obtained by extrapolation of the 
data to the thermodynamic limit [5|, |6j. This method al- 
lows for a direct investigation of the T=0 behavior in 
the whole range of parameters, circumventing thus the 
usual technical difhculties encountered in the replica for- 
malism. On the other hand, to investigate in detail the 
critical behavior |3, we formulate the Landau theory in 
the vicinity of the quantum phase transition ^ ^ . The 
consistency check of results obtained using those two ap- 
proaches allows us not only to confirm their reliability, 
but also to identify an unexpectedly narrow subregime 
near the phase boundary, in which the rapid onset of 
the glassy ordering occurs. We discuss the significance 
of our findings for the current experiments on metallic 
glasses 10]. 

We consider the random Ising model that is placed in 
a magnetic field that has both the transverse and longi- 



tudinal components, 

F = - ^ J,,S^S^-^h.S,. (1) 

<ij> i 

The random interactions Jij are chosen to be infi- 
nite range and Gaussian distributed with variance J, 
that sets the unit of energy in the model, while h = 
[h7 ,Q,h^). This model has an experimental realization 
in the LiYi_2;H0:rF4 compound that has been the subject 
of recent experiments [ill Il2| . In this insulating com- 
pound, the ground state of the magnetically active Ho 
ions is the low energy Ising doublet. In addition to that, 
the long-range nature of dipolar interactions between the 
spins enables us to perform the treatment in the large co- 
ordination limit. Disorder in the system arises from the 
fact that the substitutions of the Y atoms by the Ho ions 
are positionally random. The strong randomness leads 
to the clear observation of the spin-glass and ferro-glass 
phases at low concentration x |l3| |. 

To investigate the transition in the system described 
by the Hamiltonian J^lj we employ two methods that 
complement each other in their scope and range of ap- 
plicability. The main theoretical tool we use to obtain 
the detailed analytic behavior is based on the Landau 
theory approach[3,0]. Though attractive, this method is 
rigorously valid only close to the quantum critical point, 
so that the actual range of applicability of this approach 
is always difficult to assess. Hence, in addition, we also 
use the exact diagonalization scheme, in which one has 
to obtain the solution of H for a number of explicit real- 
izations of disorder (typically several tens of thousands). 
The procedure is implemented on finite systems of up to 
17 spins. The physical observables, such as gaps or spec- 
tral functions, are obtained along the lines of Ref. 6]. In 
this approach, no a priori assumptions are made, and 
its validity is limited by the reliability of the required 
extrapolations to the large size limit. The main reason 
for success of the previous applications of the method is 
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that for high connectivity models the numerical extrapo- 
lation to the thermodynamic limit is rather well behaved. 
Nevertheless, as we shall see and discuss later on, in the 
present study we find a certain range of parameters, where 
the previous statement does not hold. Remarkably, this 
circumstance allows us to gain new insight into the prob- 
lem. 

It is useful to characterize the parameter space by 
and h'^ , the longitudinal and transverse components of 
external magnetic field h respectively. The pure trans- 
verse field case was the subject of previous investigations 
At T=0, the existence of the quantum phase transi- 
tion was established for a value of h'^ ^ 0{J). At this 
point the s pin- spin dynamical local suceptibility becomes 
gaplessd^ll^. When the longitudinal field is turned on, 
the net longitudinal magnetization is immediately gener- 
ated and the critical point extends into a quantum crit- 
ical line h'^{h^). This line separates the two phases, in 
which the transverse degrees of freedom of spins are ei- 
ther disordered (large and h^) or spin-glass ordered 
(small h'^ and h^). As we shall show, the excitation gap 
closes at this critical line, becoming very small in some 
crossover region on the disordered side of the line. 

The Landau functional is constructed using the cumu- 
lant expansion about the quantum critical point at zero 
longitudinal field. Both the term with random interac- 
tions and the part with longitudinal field in the Hamilto- 
nian are treated as perturbations. This procedure im- 
plies that the longitudinal magnetic field is small com- 
pared to the primary microscopic energy scale h7 ^ J . 
The derivation is straightforward and leads to the follow- 
ing Ginzburg-Landau action|^] 
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Here r, being some function of h?" / J , is the parameter 
that governs the transition, while u and y are taken at the 
critical point {h^/J)c ~ 0(1). It is important to retain 
the quartic term, responsible for the RSB instability 0. 
We must insert then the mean field ansatz 



D{LOn) + f3qEASuj,„o, a = b, 



into Eq.Q and vary subsequently the free energy with 
respect to D{ujn), qsA and qab- The parametrization of 
qab depends, however, on the phase under consideration. 
In the disordered paramagnetic phase (PM) we must use 
the replica-symmetric ansatz qab — qsA, while in the 
spin glass phase (SG) the solution with a broken symme- 
try should be used|3,^3- The variational procedure is 



lengthy albeit identical to that performed in the previous 
works. As a result, we obtain that the equation deter- 
mining D{ijjn) is the same in both PM and SG phases 
and reads 
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This equation must be supplemented by 



2Z?(0)gEA + II^Ia + ^ = 



in the PM phase and 

ql^ = -[DiO)Kyy 



(3) 
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in the SG phase, to comprise the full system to be solved 
sclf-consistently. Though the exact treatment of this sys- 
tem is not possible, we can obtain the leading order of 
the correct solution close to the quantum critical point. 
We consider here only the case of T = 0, so that all the 
sums over Matsubara frequencies are substituted by the 
corresponding integrals. 

We notice first that, if y — 0, the complete solution 
is easily derived to be ji] D{uJn) = -^w^ + A^. The 
gap A^, that turns to zero right at the critical point, is 
determined using the following identity 



du 
2^ 
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In Eq. JSJ A^^ is the upper frequency cutoff and Ci is some 
constant of order unity. Let's assume that for y 7^ 
the leading approximation of D(ujn) contains the same 
square root singularity as for y = 0, and analyze how 
the last two terms in Eq.© affect the solution in the 
leading approximation. Simple inspection reveals that in 
the prelast term it is sufficient to put LUn = 0, A = 
while calculating the integral over wi. This contributes 
only to the renormalization of the coefhcient u before 
q-EA, so that u^ea Ui9ea- 

The last term requires, however, the calculation of the 
integral 



X Vl^^l +t^2 + A2, 



A2 



(7) 



that is difficult to perform exactly for arbitrary uj„ and 
A^. We need, however, only the leading behavior of this 
integral provided uj„, A 1. A simple estimate yields: 



K{A,iOn) =A + Buol + HC2/A), 



(8) 
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FIG. 1: Schematic phase diagram predicted by the Landau 
theory. The dashed hnes denote crossovers while the fuU hne 
is a critical line. 



where the constants A,B,Ci and C2 are some cut-off 
dependent functions. We see that the first term in 
the above expression renormahzes the critical value Vc 
(equal to uA'^/2tt for y = 0), while the contribution from 
the second one can be simply absorbed by the appropri- 
ate rescaling of temperature T in w^. The third term 
in Eq.® leads to the renormalization of the coefRcient 
before the A-dependent part of Eq. © . 

Similarly as in Ref. , we obtain that in the PM phase 



A2- 



r -rc + uiQEA 

1*2 ln[Cu2/ {r -rc + uiQea)] ' 



(9) 



where C,ui and U2 are again some A^^ dependent func- 
tions of the order unity. As a result of solution of Eqs.© 
and one can distinguish the following regimes on a 
(r — Tc, h) plane (see Fig^). 

(I) In this regime, in which h {r — rc)^/'^, qea is the 
smallest parameter and can be treated as a perturbation. 
As a result, we obtain with the logarithmic accuracy, that 
?EA = (k/i')/4A, a ^ {(r-r,)/M2ln(l/(r-r,)}i/2. This 
equation shows that when becomes non zero, ^ea also 
becomes finite even in the PM phase due to the finite 
magnetization along the longitudinal axis. 

The expression for the gap was first obtained in 
Ref. 01 and |0| that considered the h^=Q case. To an- 
swer the question of the region of validity of the Landau 
approach, we use the exact diagonalization method to 
obtain the gap as a function of at h^=0. The results 
are shown in Fig|21 The agreement at small values of A 
demonstrates the reliability of our methods and gives an 
indication of the size of the critical region. 

(II) This region is characterized by the condi- 
tion \r — rd^^"* h. In the leading approxima- 



FIG. 2: Gap vs. transverse field at /i^ = (open circles). 
The fitting function from Eq.@ with qea—O is plotted in the 
solid line. The arrow indicates the critical field. 



{(MiK/i2)/4u2ln(l//i4/3)}i/3^ while qEA 



tion A 

{(K/iV4)v/(u2/Mi)ln(l//i4/3)}2/3. 

(Ill) This regime, in which (rc — r)"^/^ ^ h, is the clos- 
est to the T = critical boundary. The EA order param- 
eter, that crosses over to its value in the glassy phase, is 
given by gEA = [{rc - r)/ui\ + {u2A'^/ui) In [l/A^] , with 
A [KUih'^/4:{rc - r)] - [2y{rc - r)^/3ulK'^]. From this 
expression it is easily seen that A vanishes at the critical 
line given by ft. = (8y/3)[(rc — r)/uiK]-^/^. 

Finally, in the SG phase: 

D{uJn) = -VQea/'^'^ - \^n\, QEA = {Tc - r)/ui, (10) 

resulting in a gapless form of the spectral density 
Imx(i^) oc Lo. 

We would like now to discuss the nature of the 
crossover between subregimes II and III in more detail. A 
rather surprising result, one obtains from the exact diag- 
onalization method, is that in fact the freezing transition 
of the transverse degrees of freedom takes place at the 
critical boundary fine given by h^^D oc \h^ — h^f'^'^ (see 
Fig. |3Jl . This result was verified by two different criteria: 
(i) the divergence of the spin-glass susceptibility given 
by "^^([xfoc]^) = 1' ^'^d (ii) the vanishing of the excita- 
tion energy gap of the regular part of the dynamical spin 
susceptibility, that corresponds to the so called "repli- 
con" modelig. It is notable that the infinite system size 
extrapolations for these two different freezing transition 
criteria do agree well. However, these results seem para- 
doxical since the Landau theory predicts a phase transi- 
tion boundary with a different functional form, namely, 
oc — /i^p/^ (and different curvature, see Fig. 

This paradox is resolved upon further scrutiny of the 
results from the Landau theory. To this end, it is im- 
portant to note that, in the presence of the non-zero 
longitudinal field, the critical behavior of the gap is dif- 
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FIG. 3: SG-PM phase boundary obtained with exact diago- 
nalization. Filled and open circles correspond to the two dif- 
ferent criteria (i) and (ii), respectively (see text). The solid 
line corresponds to the fitting function = — 2.2(/i^)^''^. 



ferent than at — 0. It takes a much slower, lin- 
ear form A(Jr) ~ 5r {5r is the distance to the criti- 
cal line), becoming the new effective small energy scale 
that characterizes the region III. This linear regime of 
A(5r) crosses over to the regime II, at values of Tc — r w 

{KUih'^/Af^uy^\n'^{l/h'^/^). Remarkably, this is pre- 
cisely the functional form obtained for the critical line 
(and gap closure) from the numerical calculation. 

The key point is that for systems of the size that one 
can diagonalize, the physics of the small gap is masked 
by the finite size effects, affecting thus the validity of ex- 
trapolations. This is a well known limitation of exact 
diagonalization studies that occurs in systems in which 
the small energy scales emerge Interestingly, simi- 
lar discrepancies are found in the prediction of the de 
Almeida - Thouless (AT) line in the classical Sherrington- 
Kirkpatrik (SK) model (with hF = 0) in a longitudinal 
field. Numerical calculations suggest that the critical 
temperature is Tc oc \h^ — h^\^/^, instead of the cor- 
rect result Tc oc \h^ — In the classical case, this 
can be an indication that the free energies of the ordered 
and paramagnetic phases are actually very close within a 
crossover region m the T - dia gram, equivalent to 
the region III of Fig. 1. The large finite-size effects 
were found also in the recent numerical simulations of 
the classical model [l^, and are possibly relevant to the 
anomalous behavior observed in experimental studies of 
the AT line For the quantum systems, such small 
gaps may be also difficult to observe in experiments as 
well as in numerical calculations. In contrast, in regions 
I and II, the r-dependence of A assumes a form similar 
to the zero field limit (see Eq. , except that r is shifted 
by the quantity miQea)- Since in region II (dropping log- 



arithmic corrections) i^ea ~ we conclude that the 

crossover line separating regions II and III may play a 
role of an apparent critical line, below which the gap, 
although finite, may assume unobservable small values. 

This outstanding feature, which was overlooked in pre- 
vious works, may have important consequences. For in- 
stance, it may be responsible for the peculiar observation 
of the quenching of the nonlinear susceptibility at the 
quantum critical point of the LiYi_a;Hoa;F4 series 21']. 
Another example is the electron glass model that was 
recently described in Ref. and for which essentially 
identical arguments apply. In this case the dynamical 
exponent is z = 1 , and we find that the crossover energy 
scale (corresponding to the gap in the Ising case) behaves 
as A ~ Sr^ and corresponds to a crossover temperature 
separating the Fermi liquid regime (at low T) from the 
quantum critical regime (at high T). The second power 
in Sr indicates an even broader quantum critical regime 
than in the Ising case. Such an extended quantum critical 
region may result in enhanced dissipation at low temper- 
atures, a possibility which may bear relevance for the 
puzzling absence of weak localization (interference) cor- 
rections in certain two-dimensional electron gases in the 
low density regime. 
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